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ABSTRACT 

The first stars in the history of the Universe are likely to form in the dense central regions of ~ 10 5 - 
10 6 Mq cold dark matter halos at z k 10-50. The annihilation of dark matter particles in these 
environments may lead to the formation of so-called dark stars, which are predicted to be cooler, 
larger, more massive and potentially more long-lived than conventional population III stars. Here, we 
investigate the prospects of detecting high-redshift dark stars with the upcoming James Webb Space 
Telescope (JWST). We find that all dark stars with masses up to 10 3 Mq are intrinsically too faint 
to be detected by JWST at z > 6. However, by exploiting foreground galaxy clusters as gravitational 
telescopes, certain varieties of cool (T e g < 30000 K) dark stars should be within reach at redshifts up 
to z w 10. If the lifetimes of dark stars are sufficiently long, many such objects may also congregate 
inside the first galaxies. We demonstrate that this could give rise to peculiar features in the integrated 
spectra of galaxies at high redshifts, provided that dark stars make up at least ~ 1% of the total stellar 
mass in such objects. 

Subject headings: Dark ages, reionization, first stars - dark matter - galaxies: high-redshift - stars: 
Population III 



1. INTRODUCTION 

The first stars in the history of the Universe are pre- 
dicted to form insi de ~1O 5 -1O 6 M0 minihalos at red- 
shifts z 10-50 (e.g. lTegmark et al.lll997tlYoshida et al.l 
I2003D . Due to the lack of efficient coolants in the primor- 
dial gas at these early epochs, the resulting population 
III stars are beli eved to be very massive (> 100 M ; e.g. 
lAbel et allfeOOa IBromm et a"Lll2002l), hot (eff ective tem- 
perature T e ff ~ 10 5 K: e.glBromm et al. 2001) and short- 
lived (~ 2-3 Myr; lSchaererll2002f) . Non-rotating popula- 
tion III stars with masses of 50-140 M Q or M > 260 M Q 
are expected to collapse directly to black holes, whereas 
stars with masses of 140-260 M q may produce lum inous 
pair-instability supernovae fe.g. iHeger et al.i r 2002|) . The 
latter may enrich the ambient medium with heavy ele- 
ments and initiate the transition to the normal mode of 
star formation (population I and II, with a characteris- 
tic stellar mass < 1 M©) known from the low-redshift 
Universe. The highly energetic radiation emitted from 
population III stars during their lifetimes may also have 
play ed an important role i n cosmic reionization at z > 6 
(e.g. iSokasian et ail [200l iTrenti fc Stiaveffil |2009) . An 
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observational confirmation of very massive population III 
stars would be an important breakthrough in the study 
of the star formation, chemical enrichment and reioniza- 
tion history of the Universe. 

The James Webb Space Telescop43( JWST) , sched- 
uled for launch in 2014, has been designed to study 
the epoch of the first light, reionization and galaxy 
assembly, but is not expected to be able to directly 
detect individual population III stars at z > 10. 
Searches for population III stars with the JWST would 
instead focus on the pair-instabilit y supernovae pro- 
duced at the end of their lifetimes ( Wcinm ann fc Lilly! 
120051) . or on ~10 5 -10j clusters of population 

Ill stars (IBromm et al.1 120011: iScannapieco et al. 120031 : 
ITrenti et alJ 12009b Uohnson et al.l 120091 : I Johnson! 120101 ). 
Other alternatives are to look for the spectral signatures 
of population III stars forming in po ckets of unenriched 
gas w ithin high-reds h ift ga l axies (pTumlinson fc Shull 
1 20001: iSchaererl 120021 120031: iDiikstra fc Wvithd l200"7t 
I Johnson et al.ll2008| ). or their integrated co ntribution to 
the infrared extragalactic background light (jSantos et al.1 
12002" ICoorav et al.ll2009l ). 

It has recently been recognized that annihilation of 
dark matter in the form of Weakly Interacting Massive 
Particles (WIMPs; e.g. the lightest supersymmetric or 
Kaluza-Klein particles, or an extra inert Higgs boson) 
may have generated a first population of stars with prop- 
erties very different from the canonical population III 
(jSpolvar et al.l [2008D . Because the first stars are likely 
to form in the high-density central regions of minihalos, 
annihilation of dark matter into standard model parti- 
cles could serve as an additional energy source alongside 
or instead of nuclear fusion within these objects. This 
leads to the formation of so-called dark stars, which are 
predicted to be cooler, larger, more massive and po- 
tentially longer-lived than conventional population III 

10 http://www.jwst.nasa.gov 
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stars (ISpolvar et alll200l Iloccdl2008l: IFreese et al l 



Iocco et all 12008): lYoon et al.l 120081: iTaoso et al.l . 
Nataraian et al.l 120091: IFreese et al.l 120091: ISpolvar et al.1 



20091 lUmeda et all 120091: iRipamonti et al.l I2009T T Sim- 
ilar effects have been seen in studies of t he impacts of 
dark matter upon population I and II stars (iSalati fc Silkl 
19891: iFairbairn et al.l [20081: IScott et al.l 120081 l2009al 
Casanellas fc Lopesll2009D . 
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A significant population of high-redshift dark stars 
could have important consequences for the forma- 
tion of intermedia te and supermassive black holes 
(jSpolyar et al.ll200"9T) , for t he cosmic ev olution of the pair- 
instability supernova rate ([Ioccoll20"09[ ). for the X-ray ex- 
tragalactic backg round and for the reio nization history 
of the Universe ([Schleicher et al.l [2009) . Effects such 
as these can be used to indirectly constrain the prop- 
erties of dark stars, but no compelling evidence for or 
against a dark star population at high redshifts has so 
far emerged. Here, we explore a more direct approach - 
the prospects for detection of population III dark stars 
using the JWST. 

When attempting to assess the detectability of dark 
stars at high redshifts, the expected lifespan of such ob- 
jects represents a crucial aspect. In principle, dark stars 
could live indefinitely, provided that there is ample dark 
matter available to fuel them. Dark stars are powered by 
gravitationally-contracted dark matter, pulled into their 
core as infalling gas steepens the gravitational potential 
during the formation phase. Because annihilation de- 
pletes the dark matter present within a star, ordinary 
fusion processes eventually take over as the dominant 
power source if the dark matter is not replenished. At 
this point, the dark star will essentially transform into a 
conventional population III star, albeit more massive be- 
cause the increased duration of the formation phase has 
allowed it to accrete more gas. The dark matter present 
within the star during formation will last only for a few 
million years, but these reserves may later be replen- 
ished by scattering of WIMPs on nucleons, causing them 
to lose energy and become captured in the stellar core. 
This could boost t he longevity o f dark stars substa ntially 
(llocco et al 1 12001 IFreese et all 12001 l Yoon et al l 12001 
iTaoso et al.l 120081: ISpolvar et al.ll200l Iloccoll200l . The 
ongoing replenishment of the dark matter through cap- 
ture and the resultant increase in longevity rely upon a 
number of strong assumptions and approximations, and 
the feasibility of such a mechan ism is still yet to be 
prove n by detailed calculations (jSivertsson fc Gondolol 
l201Ch . 

Whether the capture process will be efficient depends 
on two factors: the scattering cross-sections of the dark 
matter particles with ordinary nucleons, and the amount 
and density of dark matter available for capture from 
the star's surroundings. The WIMP-nucleon scattering 
cross-sections c an be constrained by direct detection ex- 
periments (e.g. ISavage et all 12001 lAhmed et all |200H 
and searches for neutrinos produ ced by annihilation in 
the Sun (e.g. [Abbasi ct al. 2009). However, the ques- 
tion of the amount of dark matter available for re fueling 
is more complicated ([Sivertsson &: Gondolol l2"0101 In a 
pristine halo, WIMP annihilations and scatterings dur- 
ing the formation stage would eventually deplete orbits 
with low angular momenta and result in a cavity of re- 
duced dark matter density in the vicinity of the dark star. 



Whether infalling WIMPs could restore the balance be- 
fore the star evolves into a supernova or a black hole may 
depend on the overall structural evolution of the minihalo 
(e.g. contraction, mergers with other halos), on tidal in- 
teractions of the dark star with subhalos, gas clouds and 
possibly other population III stars within the minihalo 
itself. Should a violent event cause the dark star to ven- 
ture far from the centre of the minihalo, the dark matter 
density would very quickly become too low to sustain fur- 
ther dark matter burning. At the current time, estimates 
of the dark star lifetime in the presenc e of capture range 
from a few t imes 10 5 to 10 10 years (e.g. lYoon et~aH l2008l: 
Hoccol [2001 iSivertsson fc Gondolol I2010D . Many of the 
mechanisms mentioned above for replenishing the dark 
matter in the centre of the dark star have moreover not 
yet been explored. In this paper, we will therefore treat 
the duration of the dark star phase as a free parameter. 

Model atmospheres and evolutionary histories of dark 
stars are presented in Sect. [5] The detectability of iso- 
lated dark stars, with and without the effects of gravita- 
tional lensing by foreground galaxy clusters, is explored 
in Sect. [3J In Sect. [H we explain how high-redshift dark 
stars can be distinguished from other objects based on 
their JWST colours, and discuss the possibility of de- 
tecting the spectral signatures of dark stars in the first 
generations of galaxies. Sect. [3] summarizes our findings. 
Throughout this paper, we will assume a ACDM cosmol- 
ogy with fl A = 0.73, Q M = 0.27 and H = 72 km s" 1 
Mpc" 1 . 

A fter this paper had been submitted, another pa- 
per (|Freese et al.l 12010ft dealing with the prospects of 
detecting high-redshift dark stars with the JWST was 
posted on the arXiv preprint server. The two studies 
differ in their assumptions concerning the masses of dark 
stars. Whereas we consider objects with masses up to 
~ 10 3 Mq, Freese et al. instead focus on the detectability 
of 'supermassive dark stars' (with masses up to 10 Mq). 

2. MODELS FOR DARK STARS 
2.1. Stellar structure and evolution 

As described in the previous section, two physical pro- 
cesses exist for bringing dark matter into a star: gravi- 
tational contraction, and capture by nuclear scattering. 
Studies of the structure, formation and evolution of dark 
stars have so far employed one of two simulation tech- 
niques: either a 'formationary' or a 'hydrostatic' ap- 
proach. In general, either physical process (or both) can 
be included with either simulation technique. Most stud- 
ies to date employing the formationary approach have 
included only gravitational contraction, whereas most 
studies following the hydrostatic approach have only in- 
cluded capture by nuclear scattering. There are however 
notable examples of both hydrostatic (jlocco et al J 120081) 
and formationary ([Spolvar et al.ll2009T) studies which in- 
clude both physical processes. 

The formationary approach follows the initial 
collapse of the pre-stellar gas cloud, the resul- 
tant gravitational contraction of its dark matter 
halo, and the subsequent formation of the dark star 
(ISpolvar et al.ll2008l: IFreese et al.l l 2"00l iNataraian et all 



20091: IRipamonti et all 120091: iSpolvar et all 120091: 
Ripamonti et all 12010ft . This approach captures the 
salient points of dark star formation and early evolu- 
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TABLE 1 

Dark star models 



WIMP mass 


M (Mf7)) a 

\ \'J ; 


R (cm) a 


loE,n(o) b 


Teff (K) a 


uurn \j / 


^sa (y r ) 


(max 


(yr) 


Atmosphere 


max(^ 


1 GeV 


106 


2.4 X 


10 14 


-0.612 


5.4 x 


10 3 


> 10 10 


2.0 x 10 5 


2.0 


x 


10 7 


MARCS 


10 




371 


2.7 x 


10 14 


-0.170 


5.9 x 


10 3 


> 10 10 


1.2 x 10 6 


1.2 


x 


10 s 


MARCS 


11 




690 


1.1 X 


10 14 


0.879 


7.5 x 


10 3 


> 10 10 


2.4 x 10 7 


5.0 


x 


10 8 


MARCS 


11 




756 


3.7 x 


10 13 


1.865 


1.0 x 


10 4 


> 10 10 


2.3 x 10 8 


5.0 


x 


10 8 


TLUSTY 


9 




793 


5.7 x 


10 12 


3.511 


3.0 x 


10 4 


> 10 10 


8.7 x 10 8 


5.0 


x 


10 8 


TLUSTY 


11 




824 


5.8 x 


10 11 


5.512 


1.1 x 


10 s 


6.2 x 10 6 


4.9 x 10 9 


6.2 


x 


10 f ' 


TLUSTY 


0.5 


100 GeV 


106 


7.0 x 


10 13 


0.458 


5.8 x 


10 3 


> 10 10 


6.0 x 10 6 


5.0 


X 


10 8 


MARCS 


6 




479 


8.4 x 


10 13 


0.955 


7.8 x 


10 3 


> 10 10 


2.2 x 10 7 


5.0 


x 


10 8 


MARCS 


11 




716 


1.1 X 


10 13 


2.895 


2.3 x 


10 4 


> 10 1() 


2.9 x 10 8 


5.0 


x 


10 8 


TLUSTY 


13 




756 


2.0 x 


10 12 


4.399 


5.5 x 


10 4 


2.2 x 10 9 


1.6 x 10 9 


5.0 


X 


10 8 


TLUSTY 


3 




787 


5.8 x 


10 11 


5.492 


1.1 x 


10 s 


5.9 x 10 6 


4.5 x 10 9 


5.9 


X 


10 6 


TLUSTY 


0.5 


lOTeV 


106 


2.2 x 


10 13 


1.463 


6.0 x 


10 3 


> 10 10 


1.7 x 10 8 


5.0 


X 


10 8 


MARCS 


2 




256 


2.2 x 


10 13 


1.846 


8.0 x 


10 s 


> 10 10 


3.1 x 10 8 


5.0 


X 


10 8 


MARCS 
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327 


2.0 x 


10 13 


2.036 


1.0 x 


10 4 


> 10 1() 


2.8 x 10 8 


5.0 


X 


10 8 


TLUSTY 


5 




399 


6.6 x 


10 12 


3.085 


2.5 x 


10 4 


> 10 1() 


2.9 x 10 8 


5.0 


X 


10 s 


TLUSTY 
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479 


2.9 x 


10 12 


3.879 


3.2 x 


10 4 


> 10 10 


1.9 x 10 9 


5.0 


X 


10 8 


TLUSTY 


2 




550 


6.0 x 


10 11 


5.307 


9.5 x 


10 4 


1.3 x 10 7 


3.6 x 10 9 


1.3 


X 


10 7 


TLUSTY 


0.5 




553 


4.8 x 


10 11 


5.503 


1.1 x 


10 s 


5.1 x 10 6 


3.9 x 10 9 


5.1 


X 


10 f ' 


TLUSTY 


< 0.5 



a from lSpolvar et al.l 1 120091) 

in units of g cm -1 s -2 

c This observability limit is set by the requirement that a single dark star should be sufficiently bright in at least one JWST filter to j 
a 5<7 detection after a 100 h exposure, if a gravitational magnification of /i — 160 is assumed (sec Sect. [3t ■ 



tion, which are primarily governed by the annihilation 
of the gravitationally-contracted dark matter. The 
strategy is not optimised for dealing with the stellar 
evolution after formation, beca use it relies on either 
full hydrodynamic simulatio ns (|Nataraian et al.l 120091 : 
iRipamonti et al.l 120091 120101) or ana l ytical approxima- 
tions to them (jSpolvar et al.l 120081 I2009D . becoming 
either too numerically-demanding or reliant upon simple 
polytropes once the star has condensed. 

The hydrostatic approach assu mes some initial 
model for ei ther a main-sequen ce dScqtt et al.l [2 008; 
Yoonetal.l 120081: iTaoso et all 120081: iScott et all 
2009al: ICasanellas fc Lopesl 120091) or prc-mam se- 
quence star (llocco et all 120081: lUmeda et all 120091 : 
ICasanellas fc Lopesl 12009). This model is then run 
through a modified quasi- hyd rostatic stellar evolution 
code (e.g. IScott et all l2009bl) and evolved with the 
inclusion of energy injection by WIMP annihilation. 
Late-stage evolution is typically dominated by the dark 
matter distribution outside the star and the resultant 
rate at which WIMPs are captured by the nuclear 
scattering process. Progress has been made in including 
some early-sta ge effects, like gray itationally-contracted 
dark matter llocco et all I2008D and gas accretion 
(|Umeda et al. 2009), but the realism of such simulations 
ultimately suffers from their inability to deal with times 
before hydrostatic equilibrium is reached. For our 
purposes, the most significant finding in hydrostatic 
studies was the possibility that the lifetimes of dark stars 
might be extended. In this situation, WIMP-dominated 
models 'stall' (either temporarily or permanently) when 
they reach an equili brium configura t ion somewhere on 
the Hayashi track (|Fairbairn et all 120081 : IScott et al 
2008t llocco et all 120081: lYoon et all 120081: ITaoso et al 



20081 ). For a given stellar mass, the position i s essentially 



depen dent only on the WIMP capture rate (|Scott et all 
l2009a|) . and does not depend strongly on whether 
simulations are started from the main or pre-main 



sequence (ICasanellas fc Lopes 200§). 

The formationary approach produces more realistic 
protostellar structures, whereas the hydrostatic one al- 
lows more accurate modelling of long-term evolution. It 
is likely that the results of more sophisticated simula- 
tions, where formation and later quasi-hydrostatic evo- 
lution are treated self-consistently and detailed capture 
calculations arc also included, would resemble a superim- 
position of the hydrostatic results upon the formationary 
ones. Depending upon the time required for dark mat- 
ter captured by nuclear scattering to become a significant 
contributor to a star's energy budget, three outcomes are 
possible: 

1 . The star may sta ll directly on the ey olutionary 
paths described bv lSpolvar et al.l (|2009T l during its 
march toward the main sequence. 

2. The star may contra ct onto the ZAMS as per the 
ISpolyar et all (|2009l ) paths, and then re-inflate as 
annihilation of captured dark matter pushes it back 
up the Hayashi track. 

3. The star might travel only some of the way to the 
ZAMS, but then be turned around and partially re- 
inflate as the captured dark matter asserts itself. 

Scenario [T] would typically be associated with a very 
quick onset of annihilation from dark matter captured 
by nuclear scattering, [5] would result if capture takes a 
very long time to assert itself, and [3] is an intermediate 
scenario. The amount of time the captured population 
takes to become significant is most sensitive to the total 
capture rate, but also to the time required for equilib- 
rium to be achieved between capture and annihilation, 
and the time required for WIMPs to thermalise inside 
the star. These quantities in turn depend very strongly 
on the adopted models for the WIMP particle (mass, an- 
nihilation cross-section and scattering cross-section) and 
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the dark matter halo (density and velocity structure) . In 
particular, these dependencies are rather degenerate; the 
impact of a short equilibrium timescale can for example 
be mimicked by a denser dark matter halo. This picture 
is further complicated by the fact that accretion may also 
continue to some degree in scenarios Q] and El because it 
may not have been halted by radiative feedback when 
the star draws close to the main sequence. The more 
massive a star becomes, the more dark matter is needed 
to keep it from contracting onto the main sequence. 

Where in the HR diagram a stalled configuration oc- 
curs depends upon the total capture rate of dark mat- 
ter; different annihilation rates are required to support 
an equilibrium structure at different locations on the 
Hayashi track. For a fixed stellar mass, greater rates of 
capture are required to support a star further from the 
ZAMS, i.e. at earlier stages in its contraction. However, 
larger stellar masses require substantially higher capture 
rates, so in realistic formation scenarios this effect will 
be countered to some degree by accretion. Which effect 
dominates depends upon the actual accretion rate during 
the formation phase. 

How long the stalling phase persists is determined by 
the timescale over which the star's hydrogen remains un- 
depleted at a given position on the Hayashi track, and 
how long capture can be realistically maintained at the 
given rate. The hydrogen burning timescale is a func- 
tion of the central temperature and density of the star 
(i.e. the position on the Hayashi track), so is therefore 
also a function of the capture rate. At low capture rates, 
nuclear burning takes over relatively quickly regardless 
of how long capture continues, whereas at high capture 
rates it can in principle remain suppressed (or even en- 
tirely absent) for an indefinite period, unless capture 
drops because the WIMP halo has been depleted. 

For very low capture rates - where stars stall only 
ver y briefly near the m ain sequence (e.g. very late on 
the lSpolvar et al.ll2009l tracks) - the longest possible life- 
time of a dark star is set by the hydrogen-burning life- 
time. This is the approximate maximum time ibum that 
it would take for all the core hydrogen to be converted to 
helium, if the star were not to contract any further. For 
normal Pop III stars, this is set by the timescale in which 
hydrogen will be depleted by CNO cycle burning, as he- 
lium burning by the triple alpha process quickly creates 
sufficient C, N and O to catalyze the cycle and allow it to 
outstrip hydrogen burning by the pp-chain. In the case 
of dark stars however, depending upon how close to the 
ZAMS any stalling phase occurs, the triple alpha process 
may not become relevant until rather late. In this case, 
the hydrogen burning lifetime will be given by the rate- 
limiting step of the pp-chain (p + p — >• d + e + + u e ), such 
that 

t - 1 (1) 

Here m p is the proton mass, T c and p c are the central 
temperature and density of the star, an d r pp is the nu- 
clear reaction rate, given in ICoxl (|2000t ). For dark stars 
stalled close to the ZAMS, the true hydrogen burning life- 
time will likely be somewhat shorter, due t o the effective- 
ness o f the CNO cycle. The simulations of iSpolvar et al.l 
(2009) did not, however, take into account evolution in 
the stellar chemical composition, so estimating the time- 



varying rate of CNO-cycle hydrogen burning in the stars 
we consider is not completely straightforward. Here, we 
choose to apply only the pp-chain limit, as this is a hard 
limit in the sense that it cannot be evaded by arguments 
about the stellar mo del- dependence of any CNO abun- 
dances we might derive. 

For very high capture rates, the longest plausible life- 
time of a dark star is instead set by the time required for 
the surrounding dark matter density to drop below that 
required to sustain the star, due to self-annihilation in 
the halo. Once the density of the halo is lowered, capture 
is reduced and the star contracts and gradually moves 
onto the main sequence. The halo density required to 
power a dark star of a give n mass and lumino sity can be 
estimated using Eq. (14) in llocco et~al"1 (|2008f) . The self- 
annihilation time isA for a halo with a given dark matter 
density can then be obtai ned by e.g. inverting Eq. (6) in 
IScott fc Siverts son (200j|). The self-annihilation times 
obtained in this manner are rather approximate, relying 
on a number of assumptions we already know to be vio- 
lated in the first stars (a constant halo self-annihilation 
rate and the constant infall of WIMPs from a uniform, 
unbound dark matter halo with an assumed velocity 
structure). A more conservative approach, which is thus 
also much less dependent on the adopted halo model, 
is to take the corresponding upper limit on stellar life- 
times to be ~100tsA, which is the strategy adopted in 
this paper. The factor of 100 results from an approxi- 
mate order of magnitude uncertainty from the constant 
annihilation rate assumptio n, and the r oughly order of 
magnitude change seen bv IScott et all (|2009aD in cap- 
ture rates from WIMP halos with alternative velocity 
distributions. 

The capture and evolutionary histories of the first stars 
are clearly very dependent upon the adopted dark matter 
particle and halo models. To parameterise these uncer- 
tainties, as a simple approximation we consider models 
stalling at var i ous po sitions on the evolutionary tracks of 
ISpolvar et all (|2009f ). These models occupy similar posi- 
tions in the HR dia gram to those com puted from ZAMS 
starting models (e.g JIocco et al.l2008f) . despite the rather 
different physical interpretation of the two scenarios. We 
allow stars to stall at various positions on these tracks 
for times of up to i sta ii — 5 x 10 s yr. This upper limit 
roughly cor responds to the highest lifetimes considered 
realistic by Iocco (2009). To also take into account the 
effects of nuclear burning and halo self-annihilation, we 
limit the plausible lifetimes for any given model to values 
below 

imax = mm(t sta ll, iburn, 100 t SA ). (2) 

The parameters for these models are listed in Table DO- 
This age span allows dark stars forming as early as z = 30 
to survive until z ks 10, and those forming as late as 
z = 10 to survive until the end of reionization (z w 6). 
In 5 x 10 s yr, the minihalos hosting dark stars are likely to 

11 Due to convergence problems with the stellar atmosphere 
models for certain combinations of parameter values, the 100 GeV 
W IMP track in Table Q] has one data points less than the track 
in ISpolvar et al.l 1 12009 )7 This does not have any impact on the 
results from the present paper, since the parameter space sampling 
of the converged models is more than sufficient to reveal the over- 
all trends in magnitude and colour evolution. The 10 TeV WIMP 
track also contains one data point more than the Spolyar et al. 
track to better bracket the 8000-10000 K divide. 
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experience many mergers with other minihalos (typically 
once every ~ 10 7 vr lGreif e"talll200l . However, it is not 
clear whether this prevents or facilitates the existence of 
long-lived dark stars. Mergers may cause the dark star 
to be ejected from the center of the halo, thereby cutting 
off the supply of WIMPs, but may also channel fresh 
dark matter into the halo center. In principle, one could 
even consid er dark stars with lifetimes in excess of 5 x 
10 8 yr (e.g. iFreese et al.l l2010h . perhaps even dark stars 
surviving into the present-day era, but since the JWST is 
probably not the optimal telescope for constraining such 
models, we will not consider them any further in this 
paper. 

Here, we consider models computed with three differ- 
ent WIMP masses (lGeV, 100 Ge V and 10 TeV) in the 
'minimal capture' approximation of lSpolvar et al.l (2009) 
in which WIMP annihilation and nuclear fusion con- 
tribute equally to the total stellar luminosity; further 
details can be found in that paper. 

2.2. Model atmosphere spectra 

To compute the expected spectral energy distributions 
(SEDs) of dark s tars, we have used the MARCS stellar 
atmosphere code (jGustafsson et al. 2008) for T e g <; 8000 
K objects and the TLUSTY code ( jHubenv fc Lanzl ll995h 
for T G ff > 10000 K. Since neither code is able to han- 
dle objects in the T e g ss 8000-10000 K range, we have 
interpolated the evolutionary tracks to produce replace- 
ment points just outside this temperature region when- 
ever possible. In a few cases, it turned out to be nec- 
essary to adopt surface gravities log(q) sligh tly different 
from those given by the lSpolvar et all (|2009ft track to get 
convergence from MARCS. Neither the omission of the 
r Gff w 8000-10000 K data points, nor the log (5) devia- 
tions have any significant impact on the results of this 
paper. 

Provided that the properties of dark matter allow the 
formation of dark stars, objects of this type are ex- 
pected to be among the first stars forming in the his- 
tory of the Universe, and therefore to have extremely 
low metallicities Z, possibly at the level given by Big 
Bang nucleosy nthesis (this impl ies Z ~ 10~ 9 , mainly 
due to Li; e.g. Ilocco et all 120071 ). All models therefore 
assume primordial abundances of H and He. For com- 
putational reasons, the MARCS atmospheres assume an 
overall metallicity of Z = 2.5 x 10~ 7 (corresponding to 
[Fe/H] = —5 with the q-enhance d abundance ratios dis- 
cussed in lGustafsson et al.l 1200 8). whereas the TLUSTY 
have been computed at Z = 0. Tests indicate that this 
minor inconsistency amounts to an uncertainty in the fi- 
nal JWST magnitudes of ~ 0.01 mag, which is irrelevant 
for the present study. 

The TLUSTY models cover the restframe wavelength 
range 0.015-300 /im, whereas the MARCS model cover 
0.13-20 /im. While parts of the spectra at wavelengths 
shortward of 0.13 /im may be redshifted into the range 
of the JWST detectors (0.6-29 /im) at z > 4, the fluxes 
in this part of the spectra are too small for detection in 
the case of cool dark stars. The wavelength coverage of 
our MARCS models is therefore more than adequate for 
our needs. 

In the case of the hotter TLUSTY dark stars, cor- 
rections for foreground absorption need to be applied 
at wavelengths shortward of Lya (0.1216/im). At these 
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Fig. 1. — The predicted apparent AB magnitudes of dark stars 
at z = 1-20 in the NIRCam/F444W filter. Each solid line corre- 
sponds to a separate dark star model from Table [T] The dashed 
horizontal lines correspond to the JWST detection limits for a 10tx 
detection of a point source after 10 4 s of exposure (thick dashed) 
and for a 5<r detection of a point source after 3.6 X 10 5 s (100 
h) of exposure (thin dashed). At z = 10—20, the dark stars are 
4—14 magnitudes too faint for detection. Hence, without the mag- 
nification boost of a foreground galaxy cluster, JWST will not be 
able to detect individual population III dark stars at the redshifts 
where they formed. Long-lived dark stars surviving until the end 
of reionization (2 Rj 6) appear somewhat brighter, but are still at 
least 2 magnitudes below the detection limit. 

wavelengths, the intergalactic medium becomes increas- 
ingly nan^trairsparent at high redshifts due to HI absorp- 
tion (|Madaulll995fl , e ventually turning almost completely 
opaque at z > 6 (e.g. iFan et aLll2006f ). To simulate this, 
all TLUSTY fluxes at A < 0.1216 /im are set to zero 
whenever z > 6, but left unattenuatcd at lower redshifts. 
While this tre atment may be too crude at lower redshifts 
(|Madaulll995|) . this is not a major concern for the current 
study, which focuses on the prospects of detecting dark 
stars at z > 6. 

As dark stars evolve along the iSpolyar et al.l ((2009) 
tracks, they eventually attain temperatures similar to 
those of conventional population III stars (~ 10 5 K) and 
photoionize large volumes of gas in their surroundings. 
The resulting HII regions will contribute emission lines 
and a nebular continuum t o the overall spectra of such 
stars (Sc haererll2002l [2003), but this is not taken into 
account by our models. Because of this, the fluxes that 
we predict for the very hottest dark stars should be con- 
sidered conservative. 

To obtain JWST magnitudes, all model spectra have 
been convolved with the transmission profiles for the 
broadband filters available for the NIRCam (0.6-5 /im) 
and MIRI (5-29 /im) instruments, and calibrated using 
the AB system. This calibration, which will be used 
throughout this paper, is based on physical fluxes and 
defined so that an object with a constant flux per unit 
frequency interval of 3631 Jy has zero AB-magnitudes 
TOab hi all filters. 

The rest-frame stellar atmosphere spectra as well as 
the AB magnitudes (as a function of redshift from z = 
to z = 20 in steps of Az = 0.5) for the dark star 
models of Table [T] are available in electronic format from: 
|http:/ /www.astro.su.se/~ez| 

3. THE DETECTABILITY OF HIGH-REDSHIFT DARK 
STARS 
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Fig. 2. — The predicted apparent AB magnitudes of dark stars at z = 6 (a) and z = 10 (b), as a function of central wavelength of the 
JWST broadband filters. Each solid line corresponds to a separate dark star model from Table [T] These lines have been colour-coded 
according to the effective temperatures of the dark stars: T cft < 8000 K (red), 8000 K < T cff < 30000 K (green) and T eft > 30000 K 
(blue). The dashed horizontal lines correspond to the JWST detection limits for a 10ct detection of a point source after 10 4 s of exposure 
(thick dashed) and for a 5<r detection of a point source after 3.6 X 10 5 s (100 h) of exposure (thin dashed). The progressively brighter 
detection thresholds at central wavelengths higher than 4.4^ (log 10 A > 0.65) reflect the lower sensitivity of the MIRI instrument (central 
filter wavelengths log 10 A > 0.65) compared to NIRCam (central filter wavelengths log 10 A < 0.65). In both panels, all dark star models 
lie significantly faintward of the detection thresholds in all filters, implying that their intrinsic brightnesses are too lo w to be detected by 
JWST. However, a magnification of fi = 160 due to gravitational lensing by a foreground galaxy cluster (see Sect. 13.211 would shift all 
models upward by 5.5 magnitudes (as indicated by the vertical arrow) and allow certain varieties of dark stars into the brightness regime 
detectable by the NIRCam instrument. This is the case for some of the T e ff < 30000 K dark stars (green and red lines) at both 2 = 6 and 
z = 10. The reason why the red lines end abruptly at 1.15 /im (log 10 A = 0.06) for z = 6 and at 2.0 fim (log 10 A = 0.3) for z = 10 is that 
the short-wavelength limit (0.13 fim) of the MARCS model spectra have entered the bluer filters at these redshifts. Since this happens at 
niAB > 38, which is a brightness regime inaccessible to the JWST, this has no impact on the present study. The sharp drop in brightness 
at A < 1.5 fim (log 10 A < 0.17) along the blue and green lines at z = 10 is due to HI absorption in the foreground intergalactic medium. 



3.1. Dark stars in random fields 

In Fig.Q] we present the AB magnitudes of all dark star 
models from Table [T] as a function of redshift (z = 1-20) 
in the NIRCam F444W filter. Also included are two esti- 
mated JWST detection thresholds for point sources, in- 
dicated by dashed horizontal lines. The first one is based 
on lOcr detections for 10 4 s exposures (thick dashed) and 
the second on 5er detections after 100 h (3.6 x 10 5 s) ex- 
posures (thin dashed). The former represent the fiducial 
JWST detection limits listed on the JWST homepagO, 
whereas the latter roughly correspond to the magnitude 
limits expected for 'ultra deep field' type observations. 

It is immediately clear that all dark star models lie sig- 
nificantly below both these detection thresholds at high 
redshifts. At z = 10-20, the intrinsic luminosities of the 
dark stars convert into apparent magnitudes that are 4- 
14 magnitudes too faint. Hence, without the magnifica- 
tion boost of a foreground galaxy cluster, JWST will not 
be able to detect individual population III dark stars at 
the redshifts where they formed. Long-lived dark stars 
surviving until the end of reionization (z « 6) appear 
somewhat brighter, but are still at least 2 magnitudes 
below the detection limit. 

Some of the dark star models in Fig. [T] exhibit F444W 
magnitudes which change very little as a function of 
redshifts. This is due to the steep the spectra of the 
T G ff > 20000 K dark stars, which attain their peak 
fluxes at rest wavelengths < 0.2/im (i.e. on the short- 
wavelength side of the F444W filter for all redshifts in 
the plotted range). As the redshift is increased, intrin- 
sically brighter parts of their spectra are redshifted into 

12 http://www.jwst.nasa.gov/ 



the F444W filter, and this almost exactly compensated 
for the increased luminosity distance. 

In Fig. [5J we display AB magnitudes for all dark stars 
from Table [1] at z — 6 and z — 10 as a function of the cen- 
tral wavelengths of all broad JWST filters. The JWST 
detection limits, defined as in Fig. [TJ are also included. 
All dark star models from Table Q] lie significantly faint- 
ward of these thresholds, regardless of which filter is used. 
However, a magnification of /i = 160 due to gravitational 
lensing by a foreground galaxy cluster (see Sect. I3.2p 
would shift all models upward by s» 5.5 magnitudes (as 
indicated by the vertical arrow) and shift certain vari- 
eties of dark stars into the brightness regime detectable 
by the NIRCam instrument (but not by MIRI). This is 
the case for several of the T e g < 8000 K dark stars (red 
lines) and a couple of the 8000 < T cff < 30000 K dark 
stars (green lines) at both z = 6 and z = 10. 

In Fig. [3] we demonstrate the need for using synthetic 
stellar atmosphere spectra (as compared to pure black 
body spectra) for predictions of this type. In Fig. [3^, 
the MARCS spectra for the 106 M (T off = 5400 K) 
dark star and the TLUSTY spectra for the 756 M Q 
(T off = 10000 K) dark star from the 1 GeV WIMP track 
in Table Q] are compared to black body spectra based on 
identical temperatures and bolomctric luminosities. At 
rest frame wavelengths of A > 0.4 /im, the shape of the 
continua are very similar for the stellar atmosphere and 
black body SEDs, but the presence of breaks in the stel- 
lar atmosphere spectra (most notably the Balmer break 
at A » 0.36 /im) will introduce substantial differences 
once these are redshifted into the JWST filters. The 
stellar atmosphere SEDs also contain a large number of 
absorption lines, but these will have far smaller effect of 
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Fig. 3. — Synthetic stellar atmosphere spectra compared to black body spectra for dark stars, a) The SEDs of the 106 Mq, T e g = 5400 
K (lower SEDs) and the 756 Mq, T e g = 10000 K (upper SEDs, multiplied by 100 to avoid cluttering) dark stars predicted for a 1 GeV 
WIMP. The black lines correspond to the synthetic stellar atmosphere SEDs generated by MARCS (for the lower, T e ff = 5400 K spectra) 
and TLUSTY (for the upper T c ff = 10000 K spectra), whereas the blue lines correspond to black body spectra generated for identical 
temperatures and bolometric luminosities. The obvious differences are the lack of breaks (most notable at RJ 0.36 /im) and absorption lines 
in the black body spectra. Due to the low but non-zero metallicity ([Fe/H] = —5) of the MARCS spectrum for the 106 Mq, T e ff = 5400 
K dark star, numerous metal lines are seen shortward of 0.36 /an in the lower SED. While these lines appear to be strong due to the 
high spectral resolution of the model spectrum, they actually have have very small equivalent widths and negligible impact on the JWST 
broadband fluxes, b) The difference between the black body AB magnitudes mbb and the synthetic stellar atmosphere AB magnitudes 
mAtmos in the NIRCam F444W filter. The different lines represent the six dark star models from Table [T] for a 1 GeV WIMP. The solid 
lines correspond to the MARCS (thick line) and TLUSTY (thin lines) dark star SEDs plotted in a). There are substantial differences 
between 77Xbb an( i m Atmos f° r th e cooler (< 10000 K) dark stars, whereas the differences for hotter dark stars are below 0.5 mag. This is 
primarily because the hotter stars become progressively more black body-like at the relevant wavelengths. For instance, the 0.36 fim break 
(which gives rise to the negative mgB — niAtmos at z > 9) is far less prominent for hotter objects. 



the broadband fluxes discussed here. In Fig. [3}), we show 
the difference between the magnitudes derived from black 
body and the synthetic stellar atmosphere spectra in the 
NIRCam F444W filter. The different lines represent the 
six dark star models from Table Q] for a WIMP mass of 
1 GeV. The solid lines correspond to the MARCS (thick 
lines) and TLUSTY (thin lines) dark star SEDs plotted 
in Fig. [3K As seen, there are substantial differences for 
the cooler (< 10000 K) dark stars, whereas the differ- 
ences for hotter dark stars are below 0.5 mag. This is 
primarily because the hotter stars become progressively 
more black body-like at the relevant wavelengths. For in- 
stance, the 0.36 /im break (which makes the black body 
spectra overpredict the F444W fluxes at z > 9) is far 
less prominent in the hotter dark stars. We conclude, 
that while black body spectra may be useful for deriv- 
ing order-of-magnitude estimates of JWST fluxes for the 
hotter dark stars, detailed stellar atmosphere models are 
required to accurately predict the fluxes for cool dark 
stars at high redshifts. 

3.2. Dark stars magnified by gravitational tensing 

Fig. Q] and [2] demonstrate that all dark stars in Ta- 
ble Q] are intrinsically too faint at z > 6 to be de- 
tected by JWST, even if extremely long exposure times 
(iexp = 3.6 x 10 5 s, i.e. 100 h) are considered. The 
only hope of detecting isolated dark stars with JWST 
at these redshifts would then be to exploit the gravita- 
tional lensing provided by a foreground galaxy cluster. 
Galaxy clusters at z » 0.1-0.6 can in principle boost 
the fluxes of high- redshift objects by up to factors ~ 100 
(e.g. lBradac et al.ll2009l : iMaizv et al.ll2010D . As shown in 
Fig. [5J this would be sufficient to lift some of the cooler 
(T eff < 30000 K) dark star models above the JWST de- 
tection threshold. For each separate dark star model, the 



entries in the max(z b s ) column in Table [T] indicate the 
maximum redshifts at which 5a detections are possible 
in at least one JWST filter after 3.6 x 10 5 s (100 h) expo- 
sures, assuming a magnification of fi — 160 (see below). 
No dark stars are detectable at z > 13, even when this 
boost due to lensing is taken into account. 

How many dark stars at z » 10 would one then expect 
to detect in a survey of a single lensing cluster? This 
depends on the magnification properties of the cluster, 
on the cosmic star formation history of dark stars and on 
their typical lifetimes r (which we treat as a free param- 
eter bounded from above by £ max , as discussed in section 
| 0<r< t max ). While gravitational lensing boosts the 
fluxes of background objects, their surface number densi- 
ties are at the same time diluted by a factor equal to the 
magnification /i. In a region with angular area 9 2 and 
magnification fi, one can show that the expected number 
of dark stars Ads m the redshift interval [z m i n , z max ] is 
given by: 

t(z) ~ r SFR(t)d os (z) 2 (l + z) 3 dt 



t{z) 



/i(z)M DS 



dz 



Nv S = c9 2 I I v .°A„ f — dtdz 

where SFR(t) is the star formation rate (in units of 
Mq Mpc~ 3 yr _1 ) of dark stars at cosmic epoch t(z), Mds 
is the dark star mass (assumed to be the same for all such 
objects) and d os (z) is the angular size distance between 
the observer and source at redshift z. In the case of flat 
ACDM cosmologies, the derivative 4^ in eq.© is given 
by: 



dt 

dz 



h (i + z) [n M (i + z f + n A } 



1/2 



(4) 



When exploring the prospects of detecting isolated 
dark stars in the high-magnification regions of a 
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foreground galaxy cluster, we have adopted MACS 
J0717.5+3745 at z = 0.546 as our lensing cluster. This 
object has the largest angular Einstein radius detected so 
far, with a relatively shallow surface mass-density profile 
which boosts t he projected area corres ponding to high 
magnifications (jZitrin et al.l 120091 12010D . Both of these 
properties combine to make this the best lens currently 
known for the study of faint objects at very high red- 
shifts. For this cluster, the angular area over which the 
magnification is in the range /i = 100-300 (with an av- 
erage magnification ~p ps 160) for sources at z s = 6-20 is 
~ 0.3 arcmin 2 . 

For the cosmic star formation history SFR(t) of dark 
stars, we have explored three different scenarios from 
iTrenti fc Stiaveliil jUg ) for single population III stars 
forming through H2 cooling in minihalos. These three 
scenarios differ in their assumptions concerning the 
amount of Lyman- Werner (hereafter LW) feedback ex- 
pected during the relevant epochs. LW radiation is emit- 
ted by hot, high-mass stars, destroys H2 and inhibits 
further population III star formation along this cooling 
channel. The very massi v e, fus ion-driven population III 
stars that iSpolvar et al.l (|2009T ) envision as the descen- 
dants of dark stars will emit copious amounts of LW ra- 
diation. However, if many of the population III stars that 
form through H2 co oling (population III. 1 s tars in the no- 
tation suggested bv lGreif fc Bromm 1120061 ) go through a 
long-lived dark star phase (during which LW fluxes are 
insignificant due to the very low effective temperatures 
of these objects), this could substantially delay the onset 
of the LW feedback compared to scenarios with no dark 
stars. Hence, the amount of LW feedback at a given 
epoch is expected to depend both on the fraction /ds of 
population III. 1 stars that at some point become dark 
stars, and the typical dura tion t of this phas e. 

The three scenarios from lTrenti fc Stiaveliil (|2009() that 
we consider are hereafter referred to as standard LW, 
reduced LW and no LW, and correspond to the cosmic 
star formation histories of population III stars formed 
through H2 in minihalos depicted in their Figs. 1, 2 and 
3, respectively. The standard LW scenario, which we con- 
sider suitable for very small /ds and short r, pushes the 
bulk of population III star formation to higher redshifts 
than the other two (possibly more suitable for large /ds 
and t). Whereas all population III star formation has 
effectively ceased by z ps 10 in the standard LW model, 
significant formation of such stars continues to lower red- 
shifts in the other two scenarios. 

In Fig. SI we plot the apparent AB magnitudes of se- 
lected dark star models versus the number of dark stars 
expected within the high-magnification regions (fi — 
100-300; Jl « 160) of MACS J0717.5+3745 per unit red- 
shift interval. These plots are based on the assumption of 
/ds = 1, i-e. that all population III stars formed through 
H 2 cooling in minihalos become dark stars, and that all 
such dark stars display similar properties (in terms of 
mass, temperature, radius and lifetime). In reality, this 
may not be very realistic. The individual merger history 
of each minihalo is likely to give rise to some variation 
in the central CD M density. Some minihalos may als o 
host multiple stars (|Turk et alJl2009HStacv et al.ll2010f ). 
which could cause some or all of the population III stars 
forming in such systems to wander out of the minihalo 



centre where dark matter annihilation is the most effi- 
cient. The predictions in Fig. [4] can, however, easily be 
rescaled to other dark star fractions by shifting the curves 
downward by a factor /ds ■ 

The left column of Fig. @] displays the results for the 
690 M , T cS = 7500 K dark star from the 1 GeV WIMP 
track and the right column the corresponding results for 
the 716 M Q , T cS = 23000 K dark star from the 100 GeV 
WIMP track. Due to the different temperatures of these 
dark stars, the AB-magnitudes in the NIRCam F444W 
filter are used in the left column and the NIRCam F200W 
filter in the right column. For both of these dark star 
models, the upper limit on the lifetime set by eq. @ is 
t max = 5x 10 8 yr. For each model, we therefore consider 
lifetimes of r — 10 6 (orange lines), 10 7 (green lines), 
10 8 (blue lines) and t — 5 x 10 8 yr (purple lines). The 
three rows of panels in Fig. U correspond to the cosmic 
population III star formation histories with standard LW 
(top row) , reduced LW (middle row) and no LW feedback 
(bottom row). The different markers within each panel 
indicate the AB-magnitudes and numbers of dark stars 
at redshifts z = 10 (triangle), z = 15 (star) and z = 20 
(square). 

In general, long dark star lifetimes r imply larger num- 
bers of dark stars at detectable brightnesses. In the 
case of the 690 M Q , T cS = 7500 K dark star (right col- 
umn), considerable numbers (> 10) of dark stars are ex- 
pected to be sufficiently bright for detection (at 5a after 
3.6 x 10 5 s) in the high-magnification regions of MACS 
J0717. 5+3745 at z « 10-11, provided that their lifetimes 
are ~ 5 x 10 s yr. This holds regardless of which of the 
three feedback scenarios is adopted. For the reduced LW 
and no LW feedback models, significant dark stars are 
expected even if the lifetime is closer to r = 10 7 yr. In 
the no LW scenario, even r = 10 6 yr dark stars can 
be detected behind MACS J0717.5+3745, but only in 
modest numbers (< 10). The situation for the 716 M , 
T c ff = 23000 K dark star (right column) is similar, ex- 
cept that this model remains detectable up to a redshift 
of z fa 13. 

Dark stars with r = 5 x 10 8 yr do not show the same de- 
cline in their expected numbers when going from z = 15 
to z = 10 as dark stars with shorter lifespan do. This 
happens because r = 5 x 10 8 yr dark stars have life- 
times that exceed the cosmic age intervals between ad- 
jacent redshift bins, allowing such objects to accumulate 
at lower redshifts, even though the cosmic star formation 
rate of population III stars is declini ng at these epochs 
for all three feedback scenarios. The ITrenti fc Stiaveliil 
(2009) models do not allow us to trace the star forma- 
tion history to epochs at z < 10, even though the star 
formation rate clearly remains non-zero at z = 10 in both 
the reduced LW and no LW feedback scenarios. In fact, 
r = 5 x 10 8 yr dark stars forming at z > 10 in these 
feedback scenarios will survive in detectable numbers to 
even lower redshifts [z ~ 6), even if one artificially sets 
their formation rates to zero at z < 10. 

In summary, cool (T c g < 30000 K) and long-lived 
( T 10 7 y r ) dark stars may well be detected at z ~ 10 
in sizeable numbers within a single, ultra deep JWST 
field if one takes advantage of the magnifying power 
of a foreground galaxy cluster. In the case of high r 
(> 10 8 yr) and cosmic star formation scenarios which 
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Fig. 4. — The number of dark stars per unit redshift interval predicted within the high-magnification regions (fi = 100—300) of the galaxy 
cluster MACS J0717.5+3745, as a function of their apparent AB magnitudes in the JWST/NIRCam F444W (left column) and F200W 
(right column) filters. In the left column, the 690 Mq dark star model with T c ff = 7500 K from the 1 GcV WIMP track has been used, 
and in the right column the 716 Mq dark star with T eff = 23000 K from the 100 GeV WIMP track. Within each panel, the differently 
coloured lines correspond to dark star life times of t = 1 Myr (orange), 10 Myr (green), 100 Myr (blue) and 500 Myr (purple). Symbols 
along these lines indicate dark star redshifts of z = 10 (triangle), z = 15 (star) and z = 20 (square). The different rows correspond to the 
standard LW (top row), reduced LW (middle row) and no LW (bottom row) star formation histories for population III stars in minihalos 
IITrenti & Stiavclli 2009J). The vertical lines within each panel represent the JWST detection thresholds for IOct detections after 10 4 s 
exposures (thick dashed) and 5cr detections after 3.6 X 10 s (100 h) exposures (thin dashed). All population III stars forming through H2 
cooling in minihalos are here assumed to go through a dark star phase with identical properties, but the results can easily be generalized 
by scaling the curves downward by the dark star fraction /ds- The lines extending to the right hand borders of the plots reflect the sharply 
decreasing fluxes in the F200W filter due to foreground HI absorption in the intergalactic medium at z > 15. 
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Fig. 5. — The JWST/NIRCam rri35e — mm vs. m2oo — "1277 
colours of T c ff < 10000 K dark stars at z = 10 (red star sym- 
bols) compared to a number of potential interlopers in multiband 
surveys: star clusters or galaxies at z = 0-15 (black dots), AGN 
template spectra at z = 0-15 (yellow dots), Milky Way stars with 
T e g = 2000-50000 K and Z = 0.001 - 0.020 (blue dots) and Milky 
Way brown dwarfs with T c g = 130-2200 K (green dots). Since the 
dark stars reside a region of this colour-colour diagram that is dis- 
connected from those occupied by these other objects, it should be 
possible to identify possible dark star candidates in deep multiband 
JWST/NIRCam surveys. 

imply significant dark star formation at z < 15, several 
dark stars may be seen even if only a minor fraction 
(jbs ~ 0.01-0.1) of all population III stars forming in 
minihalos become dark stars with temperatures in the 
detectable range. 

4. DISCUSSION 

4.1. How to distinguish isolated dark stars from other 

objects 



As demonstrated in Sect. 13.21 certain varieties of z « 
10 dark stars may be sufficiently bright and numerous 
to be detected by a JWST/NIRCam survey of the high- 
magnification regions of a foreground galaxy cluster. But 
how does one identify such objects among the overwhelm- 
ing number of mundane interlopers located in front of, 
inside or beyond the lensing cluster? 

Given the many degeneracies involved in the interpre- 
tation of broadband photometry, there may well be un- 
resolvable ambiguities in some cases. However, many of 
the cooler high-redshift dark stars should stand out in 
multiband survey data because of their unusual colours. 
This is demonstrated in Fig. [5l where we plot the colour 
indices TO356 — 771444 vs. m^oo — TI277 (based on AB- 
magnitudes in the JWST/NIRCam F200W, F277W, 
F356W and F444W filters) at z = 10 for all dark stars 
from Table Q] with T cB < 10000 K. The colours of these 
models (red star symbols) are compared to the colours 
predicted for a wide range of galaxies, star clusters, active 
galactic nuclei (AGN) and Milky Way stars. The cloud 
of black dots in Fig. [5] indicate the colours of integrated 
stellar popul ations (star clust e rs an d galaxies) gener- 
ated with the IZackriss on et al.l {2001) spectral synthesis 
model. These predictions are based on instantaneous- 
burslEE Salpeter-IMF stellar populations at redshifts 

13 This is a conservative choice, since allowing for more extended 
star formation histories would only result in a more restricted 
colour coverage for these objects 



z = 0-15 with metallicities in the range Z = 0.001-0.020, 
ages ranging from 10 6 yr up to the age of the Universe 
at each redshift and a rest-frame stellar dust reddening 
of E(B — V) = 0-0.5 mag assuming the Calzctti et al. 
(2000) extinction law. Also included in Fig. [5] are the 
expected colours of foreground stars with T e g — 2000- 
50000 K and Z = 0.001 - 0.020 in the Milky Way (i.e. 
at z = 0), based on the iLeieune et al.l (|1998D compila- 
tion of synthetic stellar atmosphere spectra (blue dots), 
and the colours of Milky W ay brown dwarf s in the 130- 
2200 K range based on the iBurrows etUI (|2003L I2006T) 
models (green dots). The yellow dot s repre sent the tem- 
plate AGN spectra of lHopkins et al.l (|2007l ) for bolomet- 
ric luminosities log 10 L^ q \/Lq = 8.5-14.0, at redshifts 
z = 0-15. Since none of these potential interlopers have 
~ TO444 and rri2oo — "7277 colours that overlap with 
those of cool z m 10 dark stars, a diagnostic diagram of 
this type can be used to cull objects that are clearly not 
dark stars from multiband survey data. However, given 
that T c g < 10000 K dark stars are unlikely to attain ap- 
parent magnitudes brighter than vtlab ~ 30 at their peak 
wavelengths (even when boosted by gravitational lensing; 
see Fig. [2p) , and can realistically only be detected in one 
or two NIRCam filters, follow-up spectroscopy of the re- 
maining dark star candidates will be required to establish 
their exact nature. This will admittedly be very challeng- 
ing, but a very coarse spectrum can possibly be obtained 
for tuab ~ 30 objects with JWST/NIRSpec (assuming 
a 3.6 x 10 5 s exposure). Objects of this type may also be 
suitable t arg ets for the 42 m European Extremely Large 
Telescop(GJ(E-ELT). 



4.2. 



The spectral signature of dark stars in high-redshift 
galaxies 

The first galaxies are expected to form at z ~ 10- 
13 in CDM halos with total masses around 1 s Mp, 
IJohnson et alJ. 12008. 2009; Grcif ct al. 2008J IRicottil 



At the time of assembly, each such object is 



[201 

likely to contain a number of minih alos in which pop ula- 
tion III stars have already formed (jGreif et al.ll2008r i. If 
some of these population III stars go through a long-lived 
(r > 10 8 yr) dark star phase, several dark stars may in 
principle congregate inside the first generation of galax- 
ies and give rise to telltale signatures in their integrated 
spectra. In Fig. |6l we display the rest - frame spectrum 
predicted by the IZackrisson et all (|2001| ) population syn- 
thesis model for a 10 8 yr old, low-metallicity (Z = 0.001, 
i.e. population II), Salpeter-IMF (mass range 0.08-120 
Mq) stellar population which has formed stars at a con- 
stant rate. This population has been assigned a stellar 
mass of 10 6 Mq, which - for a 10 8 Mq halo with baryon 
fraction /bar ~ 0.7 - corresponds to ~ 10% of the bary- 
onic mass in s tars, in rough agreement with the models of 
IRicottil (|201Ct ) for a z « 10 galaxy. The predicted NIR- 
Cam magnitudes of this object lie around tuab ~ 33- 
34 which would make it sufficiently bright for detection 
with JWST if seen through a gravitational lens with 
magnification /j, > 10. Superposed on the spectrum of 
this high-redshift galaxy is the integrated contribution 
from ten T cS = 7500 K (the 690 M© models from the 
1 GeV WIMP track) dark stars (red line). These dark 

14 http: //www. eso.org/sci /facilities / celt / 
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Fig. 6.— The rest-frame spectrum of a 10 6 M Q , Z = 0.001, 
Salpetcr-IMF stellar population (stellar mass range 0.08-120 Mq) 
which has formed stars at a constant rate for 10 8 yr (black line) 
with a superimposed contribution (red line) from ten T e ff = 7500 
K, 690 Mq dark stars (from the 1 GeV WIMP track). Despite 
making up only 0.7% of the stellar mass in this galaxy, these dark 
stars give rise to an upturn in the spectrum at rest-frame wave- 
lengths longward of 0.36^tm. In a z = 10 object, this spectral 
feature should appear longward of 3.96/xm and could serve as a 
telltale signature of cool dark stars within the first galaxies. Due 
to the low but non-zero metallicity ([Fe/H] = —5) of the MARCS 
model for the 690 Mq dark stars, numerous metal absorption lines 
are seen throughout the red spectrum. While these lines appear 
to be strong due to the high spectral resolution of the model, they 
actually have very small equivalent widths and negligible impact 
on the JWST broadband fluxes. 

stars, which contribute only 0.7% of the stellar mass in 
this galaxy, give rise to a conspicuous red bump in the 
spectrum at rest-frame wavelengths longward of 0.36/im 
(this corresponds to wavelengths longer than 3.96/xm at 
z = 10). Because of this, galaxies that contain many 
cool dark stars are expected to display anomalously red 
colours. A feature like this is very difficult to produce 
through other means. For instance, the spectrum de- 
picted in Fig. |6] cannot be attributed to dust reddening, 
since no known extinction law would allow a sharp rise 
in flux at wavelengths longward of 0.36/zm to co-exist 
with a very blue continuum at shorter wavelengths. At- 
tributing the red bump as due to thermal dust emission 
is also untenable, since this would require dust radiating 
at a temperature close to that of the dark star (7500 K), 
which is higher than the sublimation temperature of all 
known types of dust. 

In Fig. [7j we display the 771356 — 171444 vs. m 2 oo — TO277 
colour evolution (black solid line) as a function of age 
for the synthetic galaxy from Fig. [6] at z = 10. Here, 
the age runs from 10 6 yr (black triangle) up to the age 
of the Universe at this redshift (w 5 x 10 8 yr). Also 
indicated are the colours of two cool dark stars: the 7500 
K model from the 1 GeV WIMP track (blue star) and 
the 5800 K model from the 100 GeV WIMP track (red 
star). Due to their low temperatures, these two dark 
stars are far redder than the model galaxy in both colours 
plotted. As shown in Fig. [5j this is generally the case for 
T c ff < 10000 K dark stars observed in these filters. The 
galaxy is here assumed to go through a short burst of star 
formation (forming stars at a constant rate for 10 8 yr), 
after which it evolves passively. This gives a conservative 
estimate of the colour difference between galaxies and 
dark stars. Allowing a more extended star formation 



4.5 
4 
3.5 



3 




Fig. 7. — The JWST/NIRCam m 356 - m 444 vs. m 2 oo - "1277 
colour evolution as a function of age for a z = 10, low-metallicity 
(Z = 0.001), Salpeter-IMF galaxy experiencing a short burst of star 
formation (10 s yr) and passive evolution thereafter (black line). 
The black triangle indicates an age of 10 6 yr and the black square 
5 X 10 s yr (roughly the age of the Universe at this redshift). The 
star symbols indicate the colours of two cool, z = 10 dark stars from 
Tabic [1] the 7500 K, 690 Mq model from the 1 GeV WIMP track 
(blue star) and the 5800 K, 106M Q model from the 100 GeV WIMP 
track (red star). Both of these (as is the case for all T ef j < 10000 
K dark stars observed in these niters; see Fig. [5} are considerably 
redder than the colours expected for galaxies, regardless of their 
age. Dashed lines indicate how the colours of the model galaxy 
(at an assumed age of 10 s yr) would shift if it were to contain 
dark stars of either of the two types. The filled circles along the 
dashed tracks indicate mixtures at which these dark stars make 
up 1% of the stellar mass in the model galaxy. These points are 
also significantly redder than the reddest point along the galaxy 
track, indicating that a ~ 1% stellar mass fraction in dark stars 
within 2 rj 10 galaxies would be detectable through multiband 
photometry. 

episode or a star formation rate that increases over time 
would only increase the discrepancy between the colours 
of dark stars and the model galaxy. The dashed lines 
indicate how the colours of a 10 8 yr galaxy would shift 
in this diagram, if it were to harbour dark stars of the 
type considered. The filled circles along each such mixing 
track indicate the position at which the dark stars make 
up 1% of the total stellar mass. Since these points are 
significantly redder in the 771356 — TO444 colour than the 
reddest point along the standard galaxy track, a ~ 1% 
stellar mass fraction in dark stars would result in very 
peculiar colours for z = 10 galaxies and should allow such 
objects to be identified as candidate 'dark star galaxies' 
in JWST multiband survey data. 

5. SUMMARY 

In this paper, we have investigated the prospects of de- 
tecting very massive, high-redshift dark stars using the 
JWST. While individual dark stars at z > 6 will be in- 
trinsically too faint for detection, we demonstrate that 
the magnification provided by a foreground galaxy clus- 
ter will make certain varieties of long-lived (r > 10 7 yr) 
and cool (T e s < 30000 K) detectable at redshifts up to 
zwlQ. We argue that it should be possible to identify at 
least some of these dark stars in photometric NIRCam 
surveys due to their peculiar colours. If the lifetimes 
of dark stars are sufficiently long, they may also congre- 
gate during the hierarchical assembly of the first galaxies. 
We find that this could give rise to distinct signatures in 
the integrated NIRCam colours of high-redshift galaxies, 
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provided that dark stars make up at least ~ 1% of the 
overall stellar mass in these objects. 
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